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1 Introduction 



The Brandenberger-Vafa mechanism pQ is one of the few proposals within string 
cosmology for a mechanism that yields dynamical cosmological growth of three spatial 
dimensions. The idea that the most basic property of our universe could follow from 
the dimensionality of fundamental strings is very appealing. For reviews see [2J [3j H] . 

A closer examination of string gas dynamics, however, reveals certain obstacles. 
To set the stage for the present work, it was found in [5] that - assuming validity 
of the lowest-order dilaton gravity equations of motion - the dilaton can quickly roll 
to weak coupling, leaving the universe trapped in the Hagedorn phase. Moreover it 
was found that - assuming the string gas could be treated as homogeneous - there 
was no sharp dimension dependence in the string annihilation rate. So given these 
assumptions, the dynamics did not single out three dimensions as special. 

By relaxing the assumption of homogeneity, it was shown in [6] that the mechanism 
can indeed operate if one takes into account the fact that when winding strings are 
dilute enough (their mean separation is large compared to their characteristic quan- 
tum thickness) they behave semi-classically and their annihilation rates are highly 
suppressed in more than three large spatial dimensions. The conclusion of [6J was 
that if the string gas fluctuates out of the Hagedorn regime to a radiation regime in 
d effective spatial dimensions, then any remaining winding modes are indeed dilute 
enough that they generically freeze out for d > 3 and annihilate for d = 3[j] Therefore 
three dimensions is the maximum number of spatial dimensions that can grow large 
cosmologically due to an initial thermal fluctuation. 

In the present work we would like to examine whether three dimensions is also 
the minimum number of large dimensions that can result from a thermal fluctuation. 
A priori it seems easier for one or two spatial dimensions to decompactify via the BV 
mechanism, since it seems more likely that winding modes will find each other in this 
lower dimensional subspace. If the BV mechanism is indeed the reason the universe 
has three cosmologically large dimensions, we should address the question of why we 
don't observe an effectively lower-dimensional universe. 

Answers to this question often turn to anthropic arguments. However in the orig- 

1 This assumes the thermal fluctuation takes place while the string coupling is still large enough 
for string interactions to be effective. It also raises a serious issue, noted but not addressed in [5] or 
the present work, of whether a dilute string gas can evolve to give a homogeneous universe. 
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inal proposal of Brandenberger and Vafa, the authors argued that successive thermal 
fluctuations in the Hagedorn phase would eventually cause the maximum number of 
dimensions to decompactify. Here we would like to examine this possibility in more 
detail. We postulate an initial fluctuation that causes one or two dimensions to grow 
large, and investigate the likelihood of subsequent fluctuations causing a total of three 
dimensions to decompactify Thus, in contrast to our previous work [5J, we allow for 
multiple thermal fluctuations. We also relax the assumption of isotropy made in [6J . 
Besides addressing the question of how many dimensions decompactify, our results 
will shed light on rapid thermal fluctuations as a possible mechanism for overcoming 
the obstacle to the BV mechanism pointed out in [5], namely the dilaton rolling to 
weak coupling. 

Essential to our discussion is the fact that, while there might be a window of 
opportunity for three dimensions to decompactify once one or two dimensions have 
grown large, there can be no such window for more than three dimensions. That 
is, anisotropic expansion of one or two dimensions will never favor the eventual de- 
compactification of more than three dimensions. To see this, consider the case of 
d > 3 isotropic large dimensions, where the equilibrium state is radiation in d dimen- 
sions and the winding modes want to decay. The decay rate of these winding modes 
depends on the size R of the large dimensions in two ways. First, there is an en- 
hancement ~ R 2 reflecting the fact that longer strings are more likely to annihilate]^] 
Second, with W the winding number and Ax ~ \Aog R the characteristic quantum 
thickness of a string of length R, the amplitude has an impact parameter suppression 
~ exp(— ( w b Ax ) 2 ) in the d — 3 large directions transverse to string collisions. This 
exponential suppression was put forward in jH] as essential to the BV mechanism, 
in the sense that winding modes will freeze out if d > 3. Now imagine we let some 
number of these dimensions shrink while others grow. We do this to effect anisotropy, 
however we preserve the total volume (and energy) so that the equilibrium state does 
not change. This will only further (and exponentially) suppress the annihilation of 
winding modes in the smaller dimensions by increasing the impact parameters along 
the large transverse dimensions. In contrast, the annihilation amplitude in the larger 
dimensions will grow only if there are d < 3 of them and the smaller dimensions are of 
order string scale with the winding modes (as was implicitly assumed in [6j) rapidly 
colliding along them. We conclude that only d < 3 dimensions can decompactify as 
a result of anisotropic fluctuations. 

Hence studying anisotropic fluctuations within a three-dimensional subspace ac- 
2 In the T-dual picture these strings carry more momentum and are more likely to interact. 
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commodates the relevant cases for decompactiflcation. Given an initial thermal fluc- 
tuation that causes some number of dimensions to grow, our goal is to study the 
likelihood of subsequent fluctuations causing three dimensions to decompactify. To 
keep the investigation tractable we will only consider one degree of anisotropy, namely 
between one large and two small dimensions, or between one small and two large di- 
mensions. 

An outline of this paper is as follows. In section [2] we set up the dynamical 
equations and discuss the equilibrium phases of the system. In section [3] we give our 
procedure for choosing initial conditions and sampling thermal fluctuations. In section 
[4] we present our numerical results, and in section [5] we discuss their implications for 
the Brandenberger-Vafa mechanism. 



2 Dynamics 

The general setup will follow the lines of [6], which the reader may refer to for more 
details. The difference is that here we are considering two (logarithmic) scale factors, 
v and A, respectively for the d\ < 3 dimensions initially unwound and growing, and for 
the d.2 = 3 — di dimensions subsequently expanding after a fluctuation, but wrapped 
with winding modes. 

We work in type IIA string theory, on a flat toroidal background with metric (in 
a' = 1 units) 



All other dimensions are frozen at the self-dual radius. We also have the homogeneous 
shifted dilaton field ip(t) = 2<fi(t) — div{t) — d 2 X(t). When the metric and dilaton are 
coupled to matter, the equations of motion are 





i=i 



i=dl + l 



<P = ^((f 2 + dxu 2 + d 2 X 2 ) 




(2) 



Sir 2 
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Here P u , P x are the pressures along the di, d 2 dimensions respectively. Note that we 
do not consider a potential for the dilaton. There is also the Hamiltonian constraint 
(or Friedmann equation) 

E = (2tt) V^(y? 2 - d 1 i) 2 - d 2 \ 2 ) (3) 
where E is the total energy in matter. We take matter to consist of 



• Winding modes, with W denoting the winding number. Since the d\ dimen- 
sions are assumed unwound to begin with, the winding modes only wrap the 
dimensions d 2 . Winding modes evolve according to 

W = -T W (W 2 - (WO 2 ) (4) 

Here (W) denotes the equilibrium average winding and Y\y is an interaction 
rate we specify below. The contribution to the energy from winding and anti- 
winding modes is E — 2d 2 We x and their contribution to the (off-equilibrium) 
pressure is P\ = —2We x . 

• Radiation, or pure Kaluza-Klein modes, with K v and K\ denoting the mo- 
mentum numbers along the d\ and d 2 dimensions respectively. These evolve 
according to 

k v = -r Kw (Kl-{K v ) 2 ) 

k x = -T Kx (K 2 x -(K x ) 2 ) 

The energy in Kaluza-Klein modes is E = 2d\K v e~ v + 2d 2 K x e~ x and their 
contributions to the pressures are P u = 2K v e~ v ', Pa = 2K x e~ x . 

• String oscillator modes which we model as pressureless matter. 



2.1 Equilibrium phases 

In the ci-dimensional isotropic case the string gas could be in two thermodynamic 
phases. In the Hagedorn phase massive and massless modes are in thermal equilib- 
rium, the pressure vanishes and the energy in matter is conserved. At lower energy 
densities the universe is radiation dominated. This occurs when the energy density 
in d dimensions satisfies 

Pd = v 7 <PH = c d T d H +1 (6) 
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with Th the Hagedorn temperature and q the Stefan-Boltzmann constant in d di- 
mensions. This can be alternatively expressed in terms of temperatures. In a volume 
Vd with energy E, a radiation gas has temperature T d = (-Tr) 3 ^ 1 . The condition (6) 
can be then expressed as 

T d < T H (7) 

That is, the universe is radiation-dominated when the would-be radiation temperature 
falls below the Hagedorn temperature. 

In the anisotropic case, however, there is an additional possibility. Recall that we 
have three large dimensions, with di of them larger than the remaining = 3 — d%. 
Besides the Hagedorn phase, the system could be found in either a 3-dimensional or 
a <ii-dimensional radiation phase. To fix the equilibrium phase of the universe we 
generalize the condition ([7]). Given the energy of the system and the 3-dimensional 

and <ii-dimensional volumes, the lower of the two temperatures T dl = ( Cd E y d ) dl+1 and 

T3 = (^|r) 3 determines the equilibrium phase. If neither of these temperatures is 
lower than Tjj then the system is in the Hagedorn phase. For brevity, we will refer 
to these phases as TZ^, TZ 3 and H. 

When solving the equations of motion, we assume that if at any time the universe 
is found in TZ dl , then the winding modes around the remaining di dimensions are 
frozen, meaning the universe will never become effectively three dimensional. In 
essence we are assuming that an lZ dl era is thermodynamically stable and that the 
system quickly relaxes in it: all massive modes, including those propagating in the 
transverse dimensions, quickly decay to radiation in d\ dimensions. This assumption 
could be lifted through a more careful study of scale-factor fluctuations, which might 
show some probability of transitions to an equilibrium phase other than TZ dl . As 
we do not attempt such a study here, we employ the conservative assumption that a 
radiation era in any subspace is thermodynamically stable. Note that this assumption 
applies to our initial time slice, and therefore reduces the space of initial conditions 
where we could hope for three-dimensional decompactification. 

In TZ dl all energy is carried by the massless Kaluza-Klein modes K v . In thermal 
equilibrium we have 

n dl : {Kv) = 2df' (8) 
(K x ) = (W) = 

In 7Z 3 the energy is split between K v and K\, E = 2d\K v e~ v + 2d 2 K\e~~ x . Recalling 
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that d = 3 = di + c?2 we have 



: | (9) 

(K x ) = -e\ (W)=0 
o 

In the Hagedorn phase we have 

W: 12V ^L (10) 

Finally we recall the off-equilibrium interaction rates. In the radiation phases they 
are given by [HI E] 

V Kv = -e" 2 ^, T Kx = - e - 2X+ * and Y w = -e 2A+ ^. (11) 

71 7T 71 

In the Hagedorn phase these rates get multiplied by a factor 16E/9 reflecting the 
enhancement due to highly excited oscillator modes. 



3 Fluctuations, initial conditions and procedure 

We consider the following scenario: At time t — the system begins with d\ < 3 
dimensions expanding. The initial size of these dimensions is z/ > and their initial 
expansion rate is z>o > 0. All other dimensions are frozen with vanishing expansion 
rates. At time tf some number of dimensions d<i = 3 — d% fluctuate to a size X(tf) 
and begin expanding at a rate \(tf) > 0. We assume that < X(tf) < v{tf). This 
models an energy fluctuation with energy flowing from matter to gravitation. We are 
interested in the possibility of decompactification, that is, annihilation of the winding 
modes wrapping the c?2 dimensions. We will study this as a function of the initial 
conditions, the fluctuation time tf and the anisotropy parameter 

r = ?±(t f ) = e A(*/>-K*/) (12) 
R u 

Depending on the energy density, the universe could start in any of the three 
phases TZ^, 7^3, or Ti. In particular, at time t = the system could be in IZ^- 
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In that mentioned earlier, we assume that the winding modes around the 

d>2 dimensions are frozen so that no further decompactification is possible. Thus in 
order to obtain three large dimensions from successive fluctuations, at time t = the 
system should find itself in T-L. Given Uq, this puts a lower bound on the energy at 
t = 0, 

E>c dl V dl T^ +1 (13) 

At fixed volume, this essentially constraints the initial value of the dilaton as we will 
see shortly. 

One might worry that starting in H at t = leads to a further restriction when 
we require that the d\ dimensions are initially unwound. The issue is that in the 
Hagedorn phase the equilibrium winding number does not vanish, unless the ratio 
E/e u is so small that no winding is energetically allowed. In practice this is not a 
concern, as in order for the eventual decompactification of three dimensions to take 
place, the initial value of vq must be large enough, and the initial energy low enough, 
that the window for decompactification of three dimensions is always consistent with 
the condition that the di dimensions start out unwound. Another way to see this 
is to assume the contrary as follows. Suppose we start deep in the Hagedorn phase 
where the di dimensions are wound. Then when the d<i dimensions fluctuate, it is 
unlikely that they will be able to push the system out of "H and into 7^.3. Thus if the 
d\ dimensions start out wound, fluctuations can occur but will most likely just leave 
the universe trapped in the Hagedorn phase. 

At this point we review our procedure for fixing initial conditions. The maximum 
value of 1 0| consistent with the supergravity approximation is \ip\ = 1. Orienting time 
so the universe is rolling to weak coupling, we set ipo = —1. This defines our initial 
time slice; the equations of motion then guarantee that \ip(t)\ < 1 and <p(t) < for 
all t. To explore the dependence on ip we will scan over a range of values specified 
below, consistent with the universe starting in the Hagedorn phase. Given (fio and 
tfio, the expansion rate z>o is chosen at random from its Gaussian distribution 

exp^ocexpK ^^V o 2 ] (14) 

Finally we make a choice for uq, which controls the size of the dimensions d\. In 
practice we will test two values, vq = 3 and v = 5. Note that via the constraint ([3]), 
the quantities (p , (fi , z> are sufficient to determine the energy in matter at t — 0. 
Given the energy in matter and a choice for z/ , K v (t = 0) is set to its equilibrium 
value. 
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We then evolve the system to time tf using the equations of motion Q, 
At this point we effect a fluctuation for the dimensions d,2. To model a fluctuation we 
draw energy from the bath of heavy oscillator modes and re-distribute that energy to 
the other matter and metric modes, i.e. to the expansion rates. We specify the size of 
the fluctuation, that is the value of \(tf), using the anisotropy parameter r defined 



in (f2). In practice we will scan over the range 0.1 < r < 1. Then using yi(£/), y?(t/) 
and v{tf) we evaluate the energy in matter for A = and tentatively determine the 
equilibrium phase of the universe. This allows us to choose the expansion rate \{tf) 
from a Gaussian distribution e~ x2 ^ 2a2 \ where in the Hagedorn phase the variance is 

a H = 2(27r)2rf 2 ' an d m t ne radiation phase IZ3 the variance is = (^-^cp(t /) 2 j d+1 o\. 

Once \{tf) is chosen we re-calculate the energy available to matter. We still need 
to determine the values of K\(tf) and W(tf). These are chosen from a uniform 
distribution. We randomly select K\(tf) in the range (-Kseif-duai, -^equilibrium) and W(tf) 
in the range (Max{0.5, (W)}, H^if-duai)- Note that the values at the self-dual radius 
(A = 0) set a lower bound on K\ and an upper bound on W. Also the value 0.5 is 
our cutoff value for rounding down to no winding: we do not allow the fluctuation 
to result in vanishing winding, as our goal is to test whether the winding modes will 
annihilate as a consequence of interactions. 

In summary, the values <p , ip , u , K u (t = 0) and \(tf) (or r) are fixed and we 
randomly choose z> , A(t/), K\(tf) and W{tf) from the distributions given above. For 
t > tf we solve the equations of motion until the winding modes either freeze out 
or annihilate completely. In practice, when we numerically integrate the equations 
of motion, the winding modes are considered annihilated when W < 0.5 and frozen 
when Tiyiy < 0.1i>. (The use of z> here reflects the fact that the winding modes 
collide in the d\ dimensions.) This lies at the heart of our test of the BV mechanism: 
the dimensions d\ are unwound and free to expand, and we test whether the winding 
modes wrapping the remaining c?2 dimensions can stay in equilibrium as they compete 
with the expansion rate i> of the d\ larger dimensions. 

A note is in order regarding fluctuations of the scale factors. Given the Hagedorn 
equation of state S = E/Th we assume the scale factors can fluctuate randomly 
with a uniform distribution, since there is no cost in entropy. This is true at leading 
order and is central to our approach. However, the (equilibrium) entropy receives 
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corrections at large radii, given by [5J E] 



AS = lot 



l-^r(2d i )- 1 (r ?i E) 



24- 



■i e -mE 



rp-l 



2K 



(15) 



for each set of di dimensions with radius R{. These corrections are negligible when 
compared to the leading term, typically by many orders of magnitude. The probability 
distribution resulting from this correction as a function of R at constant matter energy 
is very flat in the Hagedorn phase. There is a mild dip around T]iE = 2di — 1, but 
this lies outside the Hagedorn phase except when di = 1. We ignore these corrections 
and adopt a flat distribution for scale factor fluctuations. 



4 Results 



For the scale factor uq associated with the initial fluctuation we consider two possible 
values, vq = 3 and u = 5. We will see that the results only mildly depend on the 
choice of z/ . As values of z/ larger than 5 are in a sense "too large," since they 
correspond to radii far larger than the string scale, we trust that these two values 
provide a good test of the Brandenberger-Vafa mechanism. 



Once vq is chosen, the condition (13) that the system starts in H at t = translates 



to the following condition on the value of <po'. 

<Po < -dxu - log[(2vr) dl - 2 C(il T^ +1 ] (16) 

(We're using the fact that tp = —1, while the Hubble rates in ^ are much smaller 
than 1). The above is not a severe constraint since it is consistent with the prerequisite 
of weak coupling, although it constraints the case d\ = 2 more than d\ = 1 since it 
pushes us to weaker values for the coupling. The condition that the d± dimensions 
are unwound, (W) < 0.5, translates to 

^0 > -2^q - log[9/7r] (17) 



Taken together, (16) and (17) fix the range of dilaton values we want to test. What 
happens outside this range? At strong coupling, which corresponds to less initial 
energy, the system is found in TZ^ at t = 0. As discussed above, this means additional 
dimensions will never decompactify. On the other hand at weak coupling the system 
is trapped forever in the Hagedorn phase, with all dimensions wound. 
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For the fluctuation time tf we consider three possible values, tf = 1, 10, 100 in 
string units (a' = 1). This turns out to be the most decisive parameter, with tf = 1 
giving the largest window for decompactification of three dimensions. 

To summarize the initial conditions, we fix z/q G {3,5} and tf G {1, 10, 100}. We 



then scan over < r < 1 and the range of ipo dictated by (16) and ( |17[ ). For each 
choice of initial conditions we perform 100 integrations of the equations of motion to 
sample values of z> , \(tf), K\(tf) and W{tf) from the distributions mentioned above. 

Our results are shown in Figs. [T] and [2j where we plot the number of cases for 
which the winding modes annihilate as a function of (po and r. The dependence on (po 
and r is as expected, with strong coupling (large (p ) and large fluctuations (large r) 
more likely to decompactify. The dependence on u is rather mild. The constraints 



(16), (17) push us to weaker coupling as z/ is increased, but the enhancement in the 
annihilation amplitude ~ e 2X ^ tfS) = r 2 e 2u ^^ makes up for this suppression. A striking 
feature is the dependence of the plots on tf. The reason decompactification becomes 
unlikely at large tf is the rolling of the dilaton to weak coupling. This suppresses 



the rates (11) and makes it impossible for wound strings to annihilate. Finally the 
dependence on d\ can be understood because there is more "room" for an equilibrium 
radiation phase in two dimension than in one. 

In all cases, the decisive factor for decompactification is the energy density of the 
universe. This is shown in Figs. [3] and |4| On these graphs we show contour plots 
for the number of cases decompactifying, along with contours corresponding to the 
energy density that separates ~H and TZs (thick black line) and energy density below 
which more than 90% of the cases result in three dimensions decompactifying (thin 
red line). The contours of constant decompactification probability clearly track the 
contours of constant energy density. To understand this, note that the energy density 
determines whether the universe is in the Hagedorn or radiation phase. As in [6], the 
universe may very well stay trapped in ~H, although what we mean by this here is 
that the winding modes freeze out while the d\ dimensions keep growing. In a narrow 
region - near the critical energy density separating "H and 1Z 3 - it is possible that (W) 
drops to zero while interactions are still efficient and the winding modes annihilate. 
This is responsible for the spiky ridge seen in Fig. [T] 

Since energy density is the most important parameter, in Table [T] we give the 
range of energy densities at time tf that lead to three dimensions decompactifying. 
We also express our results as a function of d\, vq and tf by giving the percentage of 
cases that decompactify three dimensions in Table [2} 
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d\=l, v =3, t f =l 




di=\, v =5, t f =l 




d\=\, V()=3, / * =10 




d\=\, V()=5, // = 10 




rf 1= l, v =3, //-!()() 




rf 1= l, v =5, // -100 




Figure 1: Number of cases decompactifying (out of 100) as a function of ipo and r, for 
Vq = 3, 5 and ty = 1, 10, 100 when cZi = 1. 



11 



di=2, v =3, t f =l di=2, v =5, t f =l 




-7 -11 



di=2, V()=3, // = 10 ^i=2, V()=5, / / = 10 




-7 -11 



Figure 2: Number of cases decompactifying (out of 100) as a function of ipo and r, for 
fo = 3, 5 and tf = 1, 10, 100 when cZi = 2. 
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dt=l, v =3, tj -I d\=\, v„=5, t f =l 




0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

r r 



d\=\, v„=3, /y-100 f/ 1= l, y =5, ^=100 




0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

r r 



Figure 3: Contours of constant probability of decompactification for d± = 1. The thick 
black line is the energy density separating the radiation and Hagedorn phases. The thin 
red line is the energy density below which more than 90% of the cases decompactify. 
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di=2, v =3, t f =l 




r 



d\=2, v =3, /y = 10 




0.2 0.4 0.6 0.8 1.0 



r 



d { =2, v„=3, //-KM) 




0.2 0.4 0.6 0.8 1.0 



d\=2, v„=5, fy=l 




0.0 0.2 0.4 0.6 0.8 1.0 



r 



d\=2, V()=5, £/=10 




d 1= 2, v„=5, // -100 




Figure 4: Same as Figure [3] but for d\ = 2. 
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di = 1 




*/ = l 


t f = 10 


*/ = 100 


v = 3 
v = 5 


(1.5 ~ 65) x 10- 9 
(2.4 ~ 640) x 10~ n 


(1.5 ~ 15) x 10~ 9 
(2.4 ~ 195) x 10- 11 


(1.5 ~ 2.8) x 10" 9 
(2.4 ~ 7.1) x 10" 11 


di = 2 




t f = l 


tf = 10 


t f = 100 


i/ = 3 
u = 5 


(1.4 ~ 6.5) x 10" 8 
(2.3 ~ 10) x 10~ 9 


(1.4 ~ 3.1) x 10~ 8 
(2.3 ~ 3.9) x 10~ 9 







Table 1: The range of energy densities at time tf that lead to three dimensional decom- 
pactification. 



di = 1 




t f = l 


t f = 10 


t f = 100 


1/0 = 3 
v Q = 5 


82 
65 


53 
47 


20 
17 


di = 2 




t f = l 


t f = 10 


t f = 100 


i/o = 3 
u = 5 


70 

43 


26 
13 


1 

0.1 



Table 2: Percentage of cases yielding three-dimensional decompactification, showing the 
dependence on d\, vq, tf. 

5 Discussion 

In this paper we studied a model in which successive thermal fluctuations are able 
to decompactify a total of three dimensions. Although these dimensions are initially 
anisotropic, it was argued in [7] that string gas cosmology will eventually lead to 
isotropization. So roughly speaking our results support the BV mechanism, as pro- 
viding a candidate stringy origin of the universe we observe. However a more refined 
statement is that in our model the probability of decompactifying three dimensions 
depends on several parameters, in particular on the fluctuation time tf and the en- 
ergy density E(tf). We now discuss the implications of these results for the BV 
mechanism. 
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We found that three large dimensions can arise provided scale factor fluctuations 
in the Hagedorn phase occur frequently, on a timescale tf ~ \/a' '. These frequent 
fluctuations are necessary so that the dilaton does not have time to roll to weak 
coupling. It seems reasonable to take tf ~ y/a', since y/a' is the natural timescale of 
the system in the Hagedorn phase. 



Assuming tf ~ \/a', the energy density E{tf) is basically fixed by the initial value 
of the dilaton. To address the dependence on this quantity consider two successive 
fluctuations, the first at time t = 0, the second at time t ~ \/a'. The initial value 
of the dilaton </?o determines the subsequent possibilities. First, there is a relatively 
narrow range towards strong coupling where the equilibrium phase at t = is IZ^- 
In that case, other than invoking anthropic arguments, we see no way to argue for a 
three-dimensional outcome. But if the initial coupling is not too strong, then we start 
in % and the following possibilities remain. One possible outcome, which fortunately 
occupies the smallest volume in the space of initial conditions, is that the system 
lands in the part of IZ3 where the winding modes wrapping the (I2 smaller dimensions 
cannot annihilate. In that case, we again expect only d\ dimensions to grow large. 
But for most initial conditions the universe either ends up in the part of TZ^ where a 
total of three dimensions decompactify, or remains stuck in the Hagedorn phase. 

This latter possibility, of staying in the Hagedorn phase, may seem discouraging. 
But in fact it could work in favor of the BV mechanism, because roughly speaking it 
brings us back to square one. As time goes by the scale factors may evolve slightly, 
but if fluctuations occur frequently enough this evolution is inconsequential. The key 
point is that even if the d\ dimensions are unwound and experience positive pressure, 
when the system is in equilibrium in % their growth over many string times is tiny. 
This is because the energy in radiation (~ \/E) is much smaller than the energy in 
matter (~ E), and furthermore the matter energy is very nearly independent of v. 
This nearly constant large energy keeps the velocity of the dilaton large, but also 
results in large "dilaton friction" and a negligible change in v. This can be seen by 
writing the equations of motion in the form 



(18) 



In our numerical solutions, even for tf = 100, the most that v changed was by one 
part in a hundred, while the energy remained constant to one part in a few thousand. 
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Thus the d\ dimensions can "hover" in the Hagedorn phase for a long time, which 
allows additional fluctuations to take place. Of course during this time the dilaton is 
rolling monotonically towards weak coupling [5], so if too much time goes by before 
another fluctuation takes place, the string coupling will become so small that the 
winding modes are unable to annihilate. But with a fluctuation timescale tf ~ \fa' ', 
this may not be a significant concern. To summarize, within our model we find that 
there is a considerable window in which scale factor fluctuations favor the eventual 
decompactification of three dimensions, provided such fluctuations occur on timescales 
~ \fa! . 

We conclude with a few issues which must be addressed in order to reach a defini- 
tive conclusion regarding the robustness of the BV mechanism. 



• We found that the BV mechanism can operate provided there are fluctuations 
into a regime where winding strings are dilute. This raises the crucial issue 
of whether a dilute string gas in the early universe can evolve towards a ho- 
mogeneous cosmology at late times. Studying this requires going beyond our 
mini-superspace approximation, in which we only kept the homogeneous modes 
of the metric and dilaton. 

• We found that the BV mechanism can operate provided fluctuations in the 
scale factors occur on timescales of order \fa' . It seems reasonable to assume 
tf ~ y/a' in the Hagedorn phase. But this assumption should be validated, for 
example by deriving the statistics of fluctuations from a study of the stochastic 
evolution of scale factors coupled to a hot string gas. 
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